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Abstract. The periodic Lorentz gas describes an ensemble of non-interacting point particles 
in a periodic array of spherical scatterers. We have recently shown that, in the limit of small 
scatterer density (Boltzmann-Grad limit), the macroscopic dynamics converges to a stochastic 
process, whose kinetic transport equation is not the linear Boltzmann equation — in contrast 
to the Lorentz gas with a disordered scatterer configuration. The present paper focuses on the 
two-dimensional set-up, and reports an explicit, elementary formula for the collision kernel 
of the transport equation. 



One of the central challenges in kinetic theory is the derivation of macroscopic evolution 
equations — describing for example the dynamics of an electron gas — from the underlying fun- 
damental microscopic laws of classical or quantum mechanics. An iconic mathematical model 
in this research area is the Lorentz gas |15j . which describes an ensemble of non-interacting 
point particles in an infinite array of spherical scatterers. In the case of a disordered scatterer 
configuration, well known results by Gallavotti |12j . Spohn [22], and Boldrighini, Bunimovich 
and Sinai [5] show that the time evolution of a macroscopic particle cloud is governed, in the 
limit of small scatterer density (Boltzmann-Grad limit), by the linear Boltzmann equation. 
We have recently proved an analogous statement for a periodic configuration of scatterers [16J, 
|17| . In this case the linear Boltzmann equation fails (cf. also Golse |13j). and the random 
flight process that emerges in the Boltzmann-Grad limit is substantially more complicated. 

In the present paper we focus on the two-dimensional case, and derive explicit formulae for 
the collision kernels of the limiting random flight process. These include information not only 
of the velocity before and after the collision (as in the case of the linear Boltzmann equation), 
but also the path length until the next hit and the velocity thereafter. Our formulae thus 
generalize those for the limiting distributions of the free path length found by Dahlqvist 
[10J, Boca, Gologan and Zaharescu [3], and Boca and Zaharescu @]. The higher dimensional 
case is more difficult, and we refer the reader to [16], [T7], [18] for further information. The 
asymptotic estimates for the distribution tails in [18J improve the bounds by Bourgain, Golse 
and Wennberg [6]. 

Our results on the Boltzmann-Grad limit complement classical studies in ergodic theory, 
where the scatterer size remains fixed. Bunimovich and Sinai [7J showed that the dynamics of 
the two-dimensional periodic Lorentz gas is diffusive in the limit of large times, and satisfies a 
central limit theorem. They assumed that the periodic Lorentz gas has a finite horizon, i.e., the 
scatterers are configured in such a way that the path length between consecutive collisions is 
bounded. The corresponding result for infinite horizon has recently been established by Szasz 
and Varju [2T] following initial work by Bleher [2]. See also the recent papers by Dolgopyat, 
Szasz and Varju [TT], and Melbourne and Nicol pjj], [20] for related studies of statistical 
properties of the two-dimensional periodic Lorentz gas. The case of higher dimensions is still 
open, even for models with finite horizon, cf. Chernov [9], and Balint and Toth pQ. Scaling 
limits that are intermediate between the Boltzmann-Grad and the limit of large times have 
been explored by Klages and Dellago [E~ 
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This paper is organized as follows. We begin by recalling the necessary details from |16| . 
[17j . then derive our main results — the explicit formulae for the collision kernels of the two- 
dimensional Lorentz gas in the Boltzmann-Grad limit. The final two sections study the im- 
plications for the distribution of free path lengths, and asymptotic properties of the collision 
kernels. 

To describe the main results of [16], [T7], let us fix a euclidean lattice £ cM. d , and assume 
(without loss of generality) that its fundamental cell has volume one. We denote by JC P C R 
the complement of the set B d + £ (the "billiard domain"), and T x (/Cp) = IC P x S^ _1 its unit 
tangent bundle (the "phase space"), with q € JC P the position and v £ S^ _1 the velocity of 
the particle. Here B d denotes the open ball of radius p, centered at the origin, and Sf- 1 the 
unit sphere. The dynamics of a particle in the Lorentz gas is defined as the motion with 
unit speed along straight lines, and specular reflection at the balls B d + £. (We will here 
also permit more general scattering processes.) In order to pass to the Boltzmann-Grad limit, 
it is convenient to rescale the length units in such a way that the mean free path length 
remains constant as p — > 0. This is achieved by introducing the macroscopic coordinates 
(x,v) = ( P d - 1 q,v)GT 1 ( P d - 1 lC p ). 

We take (xq,Vq) 6 T l {p d ~ l K, p ) as the initial data, and denote by t\(xq,vq) the time until 
the first collision (the free path length), v\{xq,Vq) the velocity thereafter, and analogously by 
Tk(xQ, vq) and v k (xQ, vq) the time between the (k — l)st and fcth collision and the subsequent 
velocity. The following theorem proves the existence of a joint limiting distribution for random 
initial data in the limit p — > 0. The statement is a variant of the more general Theorem 4.1 
in [Ej. 

Theorem 1. Fix n € N. For any Borel probability measure A on T 1 ^) which is abso- 
lutely continuous with respect to the Liouville measure vol K d x volgd-i, and for any bounded 

continuous function f : T 1 (M a! ) x (R>o x S d ~ l ) n —* R, we have 

(1) lim / f(x ,v ,T 1 (x ,v ),v 1 (x ,v ), . . . ,T n (x ,v ),v n (x ,v )) dA(x ,v ) 

= \ f(xo,vo,£ 1 ,v 1 ,...,£ n ,v n )p(v o ,£ 1 ,v 1 )p o0 + (v 1 ,^,v 2 ) 

JT 1 (M d )x(R >0 xSf- 1 )™ 

n 

'■"Po,/3+ (Vn-l,£n,V n )dA(x ,V ) d£ k dvolgd-i (v k ) . 

k=l 

The collision kernels p(vq, ^i, v{) and p p+ {v\, £2, ^2) are in particular independent of the 

choice of A and C, and are characterized in [TB], [T7] by certain measures on the homogeneous 
spaces ASL(d, Z)\ ASL(<i, R) and SL(d, Z)\ SL(<i, R), respectively. In the present paper we will 
show how, in dimension d = 2, this abstract description can be turned into explicit, elementary 
formulae; the same task seems substantially more difficult in dimension d > 2. 

Theorem [TJ forms the key ingredient in the proof of the existence of a stochastic process that 
governs the particle dynamics in the limit p — > 0, see [IT] for details. Specifically, a particle 
cloud with initial density /o evolves in time t to the density ft given by 

(2) f t (x,v)= f(t,x,v,£,v + )d£dvol s d-i(v + ), 
where / is the unique solution of the differential equation 

(3) [d t + v-V x - <%] f(t, x, v, £, v+) = / f(t, x, v ,0, v)p + (v, f , v+) dvol s d-i (v ) 

Jgd-l ' "0 1 

subject to the initial condition f(0,x,v,£,v + ) = fo(x,v)p(v,£,v+). Equation (|3|) may be 
viewed as a substitute for the linear Boltzmann equation, which describes the macroscopic 
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dynamics in a random, rather than periodic, configuration of scatterers |12| . [22], [5]. In- 
dependently of our studies in [16], [17] . Caglioti and Golse [8j have recently proposed an 
equation analogous to ([3]) in dimension d = 2, by assuming an independence hypothesis which 
is analogous to the statement of Theorem [T] above. 

In order to treat more general scattering processes than specular reflection, we recall the 
basic facts from [17} Section 2.2], limiting our attention to the case d = 2. For a scatterer 
centered at the lattice point m G C C M 2 , we use the boundary coordinates (v,w) G S\ x S}, 
where v is the velocity and w position so that q = m + pw. Let 

(4) S_ := {(v, w) G S\ x S\ : v ■ w < 0} 

be the set of incoming data at a given scatterer, i.e., the relative velocity and position at the 
time of collision. The corresponding outgoing data is parametrized by the set 

(5) S+ := {(v,w) G Si x S{ : v w > 0}. 
We define the scattering map by 

(6) Q:S-^S+, 0-,u>-) >->■ (v+,w+). 

In the case of the original Lorentz gas the scattering map is given by specular reflection, 

(7) 6(u, w) = (v — 2(v ■ w)w, w). 

Let @i(v,w) G S\ and ©2(1*, if) G S\ be the projection of Q(v,w) G S\ onto the first and 
second component, respectively. We assume in the following that 

(i) the scattering map G is spherically symmetric, i.e., if (v + ,w + ) = Q(v,w) then 
(v + K, w + K) = Q(vK, wK) for all K G 0(2); 

(ii) if w = —v then v + = —v; 

(iii) O : S- — > 5+ is C 1 and for each fixed v G S\ the map w 1— > Q±(v,w) is a C 1 
diffeomorphism from {w £ S} : v ■ w < 0} onto some open subset of S\. 

We will write (p(v,u) G [0,2-7r) for the angle between any two vectors v,u G M 2 \ {0}, 
measured counter-clockwise from v to u. Using the spherical symmetry and &i(v, —v) = —v 
one sees that there exists a constant < Bq < it such that for each v G S\, the image of the 
diffeomorphism w 1— > Qi(v,w) is 

(8) V v := {u G S\ : B @ < <p(v,u) < 2vr - B e }. 

Let us write (3~ : V v — > {if G S^ : i> • tu < 0} for the inverse map. Then /3~ is spherically 
symmetric in the sense that (3~ K (uK) = f3~(u)K for all v G S\, u G V v , K G 0(2), and in 
particular (3~{u) is jointly C 1 in v,u. 
We also define 

(9) £+(«) = e 2 (t>,/3-(u)) (t> € S}, tt € V w ). 

The map /3 + is spherically symmetric and jointly C 1 in u. In terms of the original scattering 
situation, the point of our notation is the following: Given any v_,v+ G S\, there exist 
w ,w+ G S\ such that 0(u_,iw_) = (i> + ,k; + ) if and only if Bq < (p(i)-,v+) < 2tt — B®, 
and in this case w_ and w + are uniquely determined, as w± = f3^_{v + ). For example, in the 
case of specular reflection ([7]) we have Bq = and 

v+ — V- 

(10) w + = w 



v., 



We denote the unit vectors (1, 0) and (0, 1) by e.\ and respectively. For v = (cos (ft, sin G 
S\ we set 

(11) K(V) : 



COS <p — 8111 ( 

sin 6 cos <b 
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Figure 1. The periodic Lorentz gas in "macroscopic" coordinates — both the 
lattice constant and the radius of each scatter tend to zero, in such a way that 
the mean free path length remains finite. 
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Figure 2. Scattering in the unit ball 



so that vK (v) = e.\. The classical impact parameter is related to the above by 

(12) 6(0) = &(«,«)=£- («#(«)) -ea, 

where $ = ip(v,u). The scattering cross section is o"(f?) = a(v,u) = We define the 

exit parameter correspondingly by 

(13) s(V) = s(u,v)=p+ K{u) (e 1 )-e 2 , 
with -d = <p(u, v). 

In the case of specular reflection ([7]), we have &(■#) = cos(-#/2), s($) = — cos(t9/2). 

We now turn to the main results of this paper. 
Theorem 2. For Vq,v,v+ € S\ and £ > 0, the collision kernel p p+ is given by 

(14) p Q (3 + (v, f , t>+) = cr(i>, v + ) $ (£, ~ s ( v , v o)) 
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where for any (>0, (—1, 1), 

(i ' T ( 1 + 
(15) *oK,i«,«) = -2 SO 



7T^ 



f * — max(|iu|,|z|) — 1 

\w+z\ 



1 



j/ 10 + Z / 

z/ w + z = 0, < 1 + [w| 
if w + z = 0, ^ > 1 + \w\ 



with 



(16) 



T(x) 



ifx < 
i/0 < rc < 1 
if 1 < x. 



It follows that $o(?> -2) is discontinuous at each point (£, u;, z) G S, where S is the curve 
S = {((1 + |z|) -1 ,-z,z) : -1 < z < l}; in fact $ maps any neighbourhood of such a point 
onto the interval [0,1]. On the other hand $o is continuous everywhere in the complementary 
region, (£,w,z) G (M >0 x (-1,1) x (-1,1)) \ S. 

We now turn to the proof of Theorem [21 which will involve computations similar to those 
in [231 Sec. 8]. 



Proof of Theorem^ Formula (|14[) follows from [17} Eq. (4.1) and Rem. 2.5]. By [161 Theorem 
4.4], we have 

(17) $o(£,w,z)=u y ({M€X 1 (y) : I?M n Rf = 0}) , 
with y := (£, z + it;). Here i?^ is the open rectangle 

(18) R[ z) = {(x u x 2 ) : < xi < £, z - 1< x 2 < z + 1}, 
and 

(19) Xi(y) := {M G X 1 : y G Z 2 M}, X x := SL(2,Z)\SL(2,R), 

cf. [HI Section 7.1]. To define the measure v y we express -X~i(y) as a disjoint union -Xi(y) = 



\_^ = \X\{kex,y), where 
(20) Xi(A;ci,y) 
with 



k- 1 




( P «\#) 



z + w 



(21) 



1 

v 1 



V G 



p(*) 



1 

k~ 2 m 1 



m G 



(22) 



-X"i(ei,y) 



The measure is now the Borel probability measure on X\(y) defined as (6/ir 2 ) times the 
standard Lebesgue measure dv on each component Xi(kei,y). 

Note that for each k > 2 and each M G X\(ke\,y) we have A; _1 y G Z 2 M n i?^ , since 
z, w; G (—1, 1). Hence we may replace -Xi(y) with X\(e\,y) in (fT71) . In view of (|20|) . 

and hence by (fTT|) . <3?o(<!;, -z) equals (6/7T 2 ) times the Lebesgue measure of the set of those 
v G (0, 1) for which the lattice spanned by 

(23) tn = (£,z + w) and a 2 = v(z + w) + C 1 ) 

has no point in . 

Now for any given t> G (0, 1) we have la\ ^ JS^*' for all £ G Z, by considering the first 
coordinate. It is a simple geometric fact that for any given £±,£2 G Z with |^| > 2 there exists 
some ^ G Z such that the point &xi + (sgn^)^ belongs to the closed triangle with vertices 

(z) 

0, a\,l\a\ + £20,2- Thus, if l\a\ + £20,2 G , then by convexity &n + (sgn^)^ lies inside 
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( z \ 

FIGURE 3. The rectangle and points in the lattice Zai + 7La<i 



(recall that and ai both lie on the boundary of R^). Therefore, if the lattice TLa\ +Za2 

( z) ( z) 

contains some point in R^ then in fact £\a\ + £2 a 2 S Rg for some l\ £ Z, 1% € { — 1, 1}. 
Considering the first coordinate we also get < i\ + i2V < 1, and thus (^1,^2) = (0,1) 
or (1,-1). Hence &o{£,,w,z) equals (6/-7T 2 ) times the Lebesgue measure of the set of those 
v £ (0, 1) for which both 0,2 and a\ — 0,2 lie outside R^' . 

(z) 

Considering the second coordinate we see that 0,2 and a\ — 0,2 lie outside R^ if and only if 
both v(z + w) and (1 — v)(z + w) — ^ 1 lie outside the interval (z — 1, z + 1). If z + u; > 

then this holds if and only if v(z + w) + > z + 1 and (1 — v)(z + w) — < z — 1, or 
in other words, v(z + to) > 1 + max(,z,«;) — We thus obtain the formula (|15p . The case 
z + w < is analogous. □ 

Theorem 3. For v,v + € Sj and £ > 0, i/ie collision kernel p is given by 

(24) p(«, £, «+) = cr(tf, «+) $(£, 6(«, «+)) 

w/iere 

00 /•! 



□ 



(25) <s>(£ iW )=J J $ (rj,w,z)dzdr]. 

Proof. The relation ((Ml) follows from [T71 Eq. (2.27)], and horn [T7J Remark 6.2]. 



The inner integral in (|25p evaluates to 

' 12 
12 



(26) 



$o(£,w,z) dz 



^(r 1 - 1) 

+4(r 1 -i + H)iog^^ ± ^ 



6 /t-l 



i-H)iog( ^i:g H ; 



^ (r i_l + H)(1 + lo g ( HJH)) 



+^(r 1 -i-ki)iog( T 



2H 



(l+H < ^ < 1H ' 

(e>Tqn)- 
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Figure 4. Graph of w) as a function of £ > 0, with u; = 1, 0.9, 0.8, . . . , 
(top to bottom). 

To state the result of the second integration in (|25|) . we define the auxiliary function \P(a) for 
a > as the solution to *'(a) = (a -1 - 1) logla" 1 - l| with ^(1) = 0; thus 



(27) tt(a) 



Li 2 (a) + (l-a)log(a- 1 -l) + loga-i(loga) 2 + 2 f a< 1 



L^a" 1 ) - (a - 1) log(l - a' 1 ) + log a - \ 



a > 1, 



with the dilogarithm Li2(z) = Y1T=1 W (I 2 1 — ■"■)■ ^ or £ m * ne ran S e 3 < C < I+jwl we se * 

1 — |w| \ \ 



(28) F(£,w) = *(e(l + H)) + *(£(1 - M)) - 21og£ + 2(1 + |tc7| log(^ + 
and for £ in the range = .] < £ < , i we set 

s ° l+|ui| s 1 — \w\ 



(29) io) = *(£(1 + H)) -log£ + (l - H + 2|u»|log(^- 

The density <&(£, w) in (f25|) can now be written as 



2 k 



+ \w\ 



(30) 



i-^ (C<|) 

l + ^[F(e,«,)-F(i,«;)-l] (i<£<^) 

AfG^^-G^,^)] (T^<e<T^ 

«>iir)- 



We have already seen in [TBI Remark 4.5] that <&(£, iu) is continuous for all (£,,w) € M>o x 
(—1,1). The fact that the formulae on the right hand side of (|30l) agree at the points £ = 
jTjT^i can be verified directly from (|2"Tl) - (f29"|) . using the identity Li 2 ( ) + Li2( 1 j^ ) = 



> N| -| 1 1 — \w\ 

* 2 

to M>o x [—1,1] by setting 



log(™)log(^p) + The function $(f , w ) can furthermore be continuously extended 



(31) 



*(£,±i) 



(e < I) 

l + -^[*(20-log(2£)-l] (£>§). 



s 
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Figure 5. Graph of $(£,w) as a function of w E (-1, 1), with £ = 0.5, 0.6, 
0.7, . . . , 2 (top to bottom). 



We finally remark that the function $(£, 0) was computed in |23} Prop. 3] in a different 
context; set $(£,0) = yb° x ' SL2 (4c) i n th e notation of that paper. 

Dahlqvist [10], Boca, Gologan and Zaharescu [3], and Boca and Zaharescu [I] have obtained 
explicit formulae for the limiting distribution of the free path lengths in dimension two. These 
can be recovered from the above expressions via the relations (cf. [TBI Remark 4.3]) 



(32) 



1 fi 



*o(0 = \j J $ °(£> w > z ) dw dz ' 



for the free path length between consecutive collisions, 

-l 



(33) 



®(£,w) dw = 2 / $ (v)dri 



1 



for the free path length from a generic initial point inside the billiard domain, and 

•1 



(34) 



*o(0 



$o(f,0,z) dz 



for the free path length of a particle starting at a lattice point (with the scatterer removed). 

We conclude with a brief asymptotic analysis of the collision kernels, when £ — > 00 (the 
limit £ — > is trivial). We assume in the following that £ > 1. A short calculation shows that, 
if wz < 0, we have <3?o(£j w, z) = 0. Otherwise, for wz > 0, 



(35) $ (H,w,z) = -^r 



£ + min(|t/;|, \z\) — 1 \ 6 / 1 — max(it, y) 



\w\ + \z\ 



vr2 T V 2£- (u + y) J' 



where u = £(1 — \w\), y = £(1 — \z\). Thus 



3 J (l-max(«,y))e- 1 + C»(e- 2 ) {wz > 0, y, u G [0, 1)) 
(36) $ (£,™,z) = -f < ' , , . . 

tt I (otherwise), 
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uniformly with respect to w, z £ (—1, 1) as £ — > oo. A simple integration yields for (|26p 



We substitute 77 = t£ in (|25|) and apply (j37j) to obtain 



/oo /■ 1 

(38) = A ^ ( (1 _ n 2 t 2 )(^)- 2 + 0((^)- 3 )) A 

= ^(i-^) 2 r 1 + o(r 2 ) 

if n G [0, 1), and &(£,w) = otherwise. Both (i37|) and ([38]) hold uniformly with respect to 
w £ (— 1, 1) as £ — > oo. These asymptotics can of course also be obtained by expanding (|26|) 
and ([30]) . respectively. 
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